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t^j- \ We study an extension of the Seiberg-Witten theory of 5d N = 1 supersymmetric 

Yang-Mills on R 4 x S . We investigate correlation functions among loop operators. These 
are the operators analogous to the Wilson loops encircling the fifth-dimensional circle 
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and give rise to physical observables of top ologically- twisted 5a N = 1 supersymmetric 
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Yang-Mills in the O background. The correlation functions are computed by using the 
localization technique. Generating function of the correlation functions of U(l) theory 
is expressed as a statistical sum over partitions and reproduces the partition function 
of the melting crystal model with external potentials. The generating function becomes 
a t function of 1-Toda hierarchy, where the coupling constants of the loop operators 
are interpreted as time variables of 1-Toda hierarchy. The thermodynamic limit of the 
partition function of this model is studied. We solve a Riemann-Hilbert problem that 
determines the limit shape of the main diagonal slice of random plane partitions in the 
presence of external potentials, and identify a relevant complex curve and the associated 
Seiberg-Witten differential. 
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1 Introduction 

The notion of random partitions has played a central role in recent studies on supersymmetric 
gauge theories. A highlight will be the work of Nekrasov and Okounkov [1] who used random 
partitions to derive the Seiberg-Witten solution [2] of Ad M = 2 supersymmetric gauge theories. 
Random partitions emerge therein through Nekrasov's formula |3J of the instanton sum for 
the gauge theories in the Q background [lj. This statistical model has been extended [5] to 
investigate the integrability of correlation functions of single-traced chiral observables. Such 
an extension of the Seiberg-Witten theory also becomes attractive to understand Ad M = 1 
supersymmetric gauge theories by providing a powerful tool [6J. 

The so-called melting crystal model provides a similar statistical model for 5d J\f — 1 
supersymmetric gauge theories [3 13 [9] and A- model topological strings [10]. This is a statistical 
model of random plane partitions, also interpreted as a g-deformation of random partitions 
[H [H], [12] , that stems from Nekrasov's formula [3] for 5d supersymmetric gauge theories in the 
Q background. 

Recently, an extension of this melting crystal model with external potentials was proposed, 
and its integrable structure was elucidated [131 E] • It was argued therein that the external 
potentials are related to loop operators of a 5d M = 1 supersymmetric Yang-Mills theory 
(SYM). The goal of this paper is to present the detail of this result along with some implications. 
In particular, we study the thermodynamic limit of the partition function of this model, and 
solve a Riemann-Hilbert problem that determines the limit shape of the main diagonal slice of 
random plane partitions in the presence of external potentials. We can thus eventually identify 
a relevant complex curve and the associated Seiberg-Witten differential. 

This paper is organized as follows. Section 2 starts with a brief review about bd M = 1 
SYM in the Q background. We introduce loop operators Ok (k = 1, 2, • • ■ ) of this theory. 
Computation of correlation functions among these operators is presented by using the localiza- 
tion technique. Generating function of the correlation functions of U(l) theory is expressed as 
a statistical sum over partitions and reproduces the partition function of the aforementioned 
melting crystal model. For an application of the localization theorem, we use equivariant de- 
scent equations of the loop operators, which is proved in Appendix B, and a related T 2 -action 
on the gauge theory is argued in Appendix A. Section 3 is concerned with a common integrable 



structure of 5d M = 1 SYM in the Q background and melting crystal model. The loop operators 
are converted to the external potentials of the melting crystal model. This eventually shows 
that, by regarding the coupling constants t = (ti,t 2 , • • •) of the loop operators as a series of 
time variables, the generating function becomes a r function of 1-Toda hierarchy. In Section 4, 
we introduce an energy functional. This functional is a quadratic form, and obtained from the 
logarithm of the statistical weight in the partition function. The formula is stated in a general 
setting and is proved in Appendix C. By using the formula, we argue the thermodynamic limit 
of the melting crystal model or the g-deformed random partitions. In Section 5, the thermo- 
dynamic limit is reformulated as a Riemann-Hilbert problem to obtain an analytic function on 
C* — /, where C* is a cylinder and / is an interval in the real axis, and satisfies suitable con- 
ditions. In Section 6, we solve the Riemann-Hilbert problem. The relevant complex curve and 
the associated Seiberg-Witten differential are presented. The vev's of the loop operators are 
particularly expressed as residue integrals of the Seiberg-Witten differential. The solution of 
the Riemann-Hilbert problem can be further analyzed by applying the classical Jensen formula 
in complex analysis. The case of the single coupling constant t = (ti,0,0, •••) is described 
as an example. The related computation is attached to Appendix D. Section 7 is devoted to 
conclusion and discussion. 

2 Loop operators of bd M = 1 SYM in Q background 

We first consider an ordinary 5d M = 1 SYM on R 4 x 5 1 . Let E be the SJJ (./V")-bundle on R 4 
with 02(E) = n > 0. A gauge bundle of the theory is the SU(N)-bundle n*E on R 4 x S 1 pulled 
back from IR 4 , where n is the projection from R 4 x S 1 to R 4 . All the fields in the vector multiplet 
are set to be periodic along S 1 . Among them, the bosonic ingredients are a 5d gauge potential 
Am{%, t)dx M and a scalar field (p(x,t) taking the value in su(N), where the coordinates x M 
represent coordinates x = (x 1 , x 2 , x 3 , x A ) of R 4 and a periodic coordinate t of S 1 . These bosonic 
fields describe a 5d Yang-Mills-Higgs system. The gauge potential can be separated into two 
parts A^x, t)dx^ and A t (x,t)dt, respectively the components of the R 4 - and the ^-directions. 
Let Ae be the infinite dimensional affine space consisting of all the gauge potentials on E. The 
4d component A^{x ) t)dx^ 1 describes a loop A(t) in Ae, where the loop is parametrized by the 
periodic coordinate of the fifth-dimensional circle. As for A t (x, t), together with (p(x,t), the 



combination A t + icp describes a loop (pit) in fi°(R 4 ,acLE ® C), which is the space of all the 
sections of adi?(g)C, where adi? is the adjoint bundle on M 4 with fibre su(N). Taking account of 
the periodicity, the same argument is applicable to the gauginos as well. The vector multiplet 
thereby describes a loop in the configuration space of the Ad theory. By using such a loop, we 
may describe bd M = 1 SYM. In the case of the Yang-Mills-Higgs system, the loop A(t) gives 
a family of covariant differentials on E as dA(t) — d + A(t). For the loop (j)(t), since it involves 
A t (x,t), we conveniently introduce a differential operator TC(t) by 

H{t) = ± + <Kt). (2.1) 

2.1 bd J\f = 1 SYM in Q background 

Via the standard dimensional reductions, 6d J\f — 1 SYM gives lower dimensional Yang-Mills 
theories with 8 supercharges, including the above theory. Furthermore, the dimensional reduc- 
tions in the Q background provide much powerful tools to understand these theories [I]. The 
Q background is a 6d gravitational background on R 4 x T 2 described by a metric of the form: 

4 

ds 2 = /^(dx^ — 2, Vadx a ) + y, (dx a ) , where two vectors V^,Vq generate rotations on 

fj.=l a=5,6 a=5,6 

two-planes (x 1 ,^ 2 ) and (x 3 ,x 4 ) in M 4 . By letting V\ = x 2 -^ — x 1 -^ and V 2 = x 4 -^ — x 3 -£^, 
they are respectively the real part and the imaginary part of the combination 

V eue2 = € 1 V 1 + e 2 V 2 , ei,e 2 eC (2.2) 

The above combination is expressed in component as V ei>e2 = VL^ u x u -£^. 

To see the dimensional reduction in the ^-background, let us first consider the bosonic part 
of the 5d SYM. The corresponding Yang-Mills-Higgs system is modified from the previous one. 
However, the system is eventually controlled by replacing TC(t) with 

7U ltea {t)=n(t)+K eue2 (t). (2.3) 

Here /C ei)€2 (t) is an another differential operator given by [15] 

K«*(t) = KleMt), + \wj»v , (2.4) 

where J^ v denote the SO (4) Lorentz generators of the system. The above operator generates 
a T 2 -action by taking the commutators with d,A(t) an d TC(t). For instance, we have 

[d A (t), ^£1,62 (t)\ = -i^v ei ,e 2 F A(t) , (2.5) 

3 



where the right-hand side means —1 times the contraction of the curvature two-form F^ = 
dA(t) + A(t) A A(t) with the vector field V eij62 . As is argued in Appendix [A], the infinitesimal 
rotation 5x^ = —V£ e2 can generate a T 2 -action on Ae- The right-hand side of (12.51) is precisely 
the infinitesimal form ( 1A.10I) of the T 2 -action. 

The supercharges Q aa and Q2 are realized in a way different from the case of ei = £2 = 0. 
Note that we use the 4d notation such that a, a and a denote the indexes of the Lorentz 
group SU{2)l x SU{2)r and the R-symmetry SU{2)i. According to the argument [IE], we may 
interpret the 5d SYM as a topological field theory [17]. Actually, by regarding the diagonal 
SU{2) of SU{2) R x SU(2)j as a new SU(2) R , we can extract a supercharge that behaves as a 
scalar under the new Lorentz symmetry. We write the scalar supercharge as Q ei ,e 2 - The gaugino 
acquires a natural interpretation as differential forms, r](x, t), ip^x^t) and ^(x, t). These give 
fermionic loops, 7](t),ip{t) an d £(£)• The main part of the Q-transformation can be read as 

Q ei , e2 A{t) = ip(t) , Q ei ,eMt) = [ d A(t), n eue2 {t)], (2.6) 

Q ei , ea W ei , ea (t) = 0, (2.7) 

where ip(t) is a fermionic loop in fi 1 (M 4 , ad-E). 

The action of the 5d SYM can be now written in a Q-exact form as 

4n 2 n 



S- 



£1,62 



fdt+f d*{Q ei , ea ,W(*)}. (2.1 



5dSYM — o , » 

9 Js 1 JSi 

The above action turns out to be a BRST gauge fixed action of the topological term ^P- J gl dt, 
where Q ei ,t 2 plays a role of the BRST charge. In this formulation, the original vector multiplet is 
splitted to several Q-multiplets, including A(t),ip(t) as the doublet and 7Y 6l , e2 (t) as the singlet. 
We may integrate out all the multiplets except the two. The integrations impose constraints 
on the two multiplets. As for the doublet, the constraints turn to involve three equations; 
A(i) = 0' \dA(t)' l P(t) ) = 0, where the symbol (+) means self-dual part of two-forms on M 4 , 
and d* A , t *.if)(t) = 0. The first equation is the anti-self dual equation for gauge potentials on E. 
The second equation is regarded as a linearization of the first. For the singlet, the constraint 
leads to the equation 

n A{t) <f>(t) + 2m 2 + d* A{t) (iv^F m + ^) = , (2.9) 

where ip(t) 2 = ip^(t)i/j^(t). In the above, d A , t s is the formal adjoint of d,A(t) an d Da(() = 
(—l)d A , t sdA(t) is the scalar Laplacian. 
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2.2 Loop operators and their correlation functions 

Taking account of the relation <p(x,t) = A t (x,t) + i(p(x,t), the following path-ordered integral 
becomes an analogue of holonomy of the 5d gauge potential: 

W {0) {x]t u t 2 ) = Pexpf- / dt<f)(x,t)) , (2.10) 

where the symbol means path-ordered integration, more precisely, it is defined by the differential 
equation 

(-!L + < f>(x,t 1 ))WW(x;t 1 ,t 2 ) = Q, W^(x;t 2 ,t 2 ) = 1. (2.11) 

at i 

The trace of the holonomy along the circle defines a loop operator by 

O {0 \x) =TtW {0 \x;R,0), (2.12) 

where R is the circumference of S 1 . The above operator is an analogue of the Wilson loop 
operator along the circle. Unlike the case of e x = e 2 = [18], it is not Q-closed except at x — 0. 
To see this, note that the Q-transformations ( I2.6P and ( 12. 7j) imply Q eui2 (f>(t) = —iy t i>(t). By 
using this, we obtain the transformation 

Q euea O<> \x)=J dfTr{wV>(x; R,t) t Vti ^{x,t)W^{x; t,0)}. (2.13) 

Since the right-hand side of the above formula vanishes only at x — 0, this means that O^ \x) 
becomes Q-closed only at x — 0. 

The above property may be explained in terms of the equivariant de Rham theory. To see 
this, let us first generalize the path-ordered integral ( 12.101) by exponentiating the combination 
^A{t) ~ i>(t) + </>(t) in place of <f>(t) as 

W(x; t!,t 2 ) = Pexpj- J ' dt(F A(t) - i/;(t) + <j>(t))(x)} , (2.14) 

where the right-hand side is the solution of the differential equation 
dW(x;t 1 ,t 2 



+ [F A{tl) (x)-ip(x,t 1 ) + (j)(x,t 1 ))AW(x;t 1 ,t 2 ) = 0, 

W(x;t 2 ,t 2 ) = 0. (2.15) 



This particularly means that W has several components, according to degrees of differential 
forms on M 4 , as W = W^ + W^ + • • ■ + W^ 4 \ where the indexes denote the degrees. We then 
generalize the loop operator (I2.12p as 

O(x)=TrW(x;R,0). (2.16) 

In parallel with W, we have the decomposition O = O (0) + C (1) H h (4) . Explicitly, £> w 

are operators of the following form: 

{1) {x) = J dtTrlw {0 \x;R,t)i;(x,t)W {0 \x;t,0) 

<-R 



{2) (x) = I dtTi(w io \x;R,t)(-)F A{t) (x)W io \x;t,0) 



+ J dt x f ' dt 2 Ti^W w (x; R,t 1 )i;(x,t 1 )W w (x;t u t 2 )i;(x,t 2 )W (0 \x-,t 2 ,0) 
i (2.17) 

Eq.f l2.13p can be now expressed as Q^^O^ = Ly e e 0^ l \ which is actually the first equation 
among a series of the equations that O^ obey. Such equations eventually show up by expanding 
component-wise the identity 

(d ei , ea + Q eil6a )O(x) = 0, (2.18) 

where d eitE2 = d — Lv £ e is the T 2 -equivariant differential on M 4 . The above identity implies in 
components T 2 -equivariant descent equations of the form 

dO^- 1 \x)+Q eu€2 O^(x)-i Vn ^ l+1 \x)=0 } 0<*<4, (2.19) 

where O^ 1 ^ = 0^> = 0. We provide a proof of the formula ( 12.181) in Appendix IB1 

We can also consider the loop operators encircling the circle many times. Correspondingly, 
we introduce 

O k (x) = TrW{x;kR,0), k = 1, 2, • • • . (2.20) 

These satisfy 

{d ei , ea + Q euea )O k (x) = 0. (2.21) 



Let us examine the correlation functions ( Y[ a J K 4 @k a ) €1 ' £2 - Since the integral / R 4 Cfc = 
J R4 Oj. is Q-closed by virtue of the formula (j2.21| . the correlation function has an interpretation 
in the topological field theory. In particular, we may compute the correlation functions by a 
supersymmetric quantum mechanics (SQM) which is substantially equivalent to the 5d SYM 
as the topological field theory. Such a SQM turns out to be an equivariant SQM on M. n [3], 
where Ai n is the moduli space of the framed n instantons. 

We may regard M. n as a Riemannian manifold. Actually, by taking a local gauge slice of 
the anti-self dual gauge potentials, the metric is induced from Ae, where Ae is endowed with 
a gauge invariant metric of the form, G(6A, 5A) = J R4 Tt5A(x) A *5A(x). Take a loop m{t) = 
(m J (t)), where m 1 denote local coordinates of M. n . We also introduce its fermionic partner 
x{t) = (x 1 ^)). Both are parametrized by the fifth-dimensional circle. The Q-transformation 
(12. 6p yields a transformation between them. To describe the transformation, note that the 
infinitesimal variation 5 A = —iv e E Fa preserves the anti-self dual equation F^ = 0, thus it 
gives a vector field on Ai n , which we denote by V ei)t2 . It is actually a Killing vector. This 
apparently follows since the aforementioned variation preserves the gauge invariant metric G. 
The Q-transformation (12.61) is eventually converted to the transformation 

Qe U eMt) = X(t) , Qe U e 2 x(t) = ~^p + V 61iea (m(*)) . (2.22) 

The supersymmetric Yang- Mills is described effectively in terms of m(t) and x{^)i conse- 
quently reduces to a quantum mechanical system on M. n - Actually, the corresponding action 
can be obtained from (12.81) by integrating out the irrelevant fields. This yields the action 

Sl)f = ^J si dt + J s dt {q ei , 2 , \Gu{^{t))x\t) d ^^ , (2.23) 

where Gu is the Riemannian metric on M. n - The above action together with the supersymmetry 
(12.221) is familiar in the physical proof [19] of the equivariant index formula for Dirac operator. 
In particular, the partition function becomes eventually the T 2 -equivariant index for the Dirac 
operator on Ai n [3]. 

The combination F^m —ip(t) +<f>(t) in ( 12.141) can be identified with a loop space analogue of 
the T 2 -equivariant curvature T^^ of the universal connection [20], where the universal bundle 
becomes equivariant by the T 2 -action on Ae x M 4 . In the computation of the correlation 
function, by virtue of the supersymmetry ( 12.221) . only the constant modes m , Xo contribute 
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to the observables. For such constant modes, the above combination is precisely identified 
with the equivariant curvature J r ei . t2 [4J. This means that Ok(x) substantially truncates to the 
T 2 -equivariant Chern character Tr e~ e i> e i. Thus we obtain the following finite dimensional 
integral representation of the correlation functions: 

II / °*. T'" = W / ^(^t 6l , e2 , M n ) n / Tr e- k « R ^ 2 . (2.24) 

a JR* I n-tnstanton fimR) ^^ J Mn a ^ 

where A T 2(- ,Ai n ) is the T 2 -equivariant A-genus of the tangent bundle of Ai n , and t £lj62 is a 
generator of T 2 that gives the Killing vector V eij(E2 . 

Introducing coupling constants t = (t\, £2, ■ ■ ■ ), the generating function of the correlation 
functions is given by 2 ei>e2 (t) = (e^ ktk -l* 4 ° k ') ' . Since n-instanton contributes with the weight 
(RA) 2nN , where A is the dynamical scale, letting Q = (RA) 2 , we can express the generating 
function as 



Z eue2 (t) = J2Q nN (e^y° k \ eue2 . (2.25) 



n=0 



2.3 Application of localization technique 

The right-hand side of the formula (j2.24p is eventually replaced with a statistical sum over 
partitions. To see their appearance, note that the integration localizes to the fixed points of 
the T 2 -action. However, the fixed points in M. n are small instanton singularities since the 
variation 5 A = — ty e E Fa vanishes there. These singularities can be resolved by instantons on 
a non-commutative M 4 . By using such a regularization via the non-commutativity, the fixed 
points get isolated, so that they are eventually labelled by using partitions [21]. 

A partition A = (Ai, A 2 , • • • ) is a sequence of non-negative integers satisfying \ > \ i+ i for 
alH > 1. Partitions are identified with the Young diagrams in the standard manner. The size 
is defined by |A| = Xl«>i ^j> which is the total number of boxes of the diagram. 

Let us examine the formula (12.241) for the U(l) theory. The relevant computation of the 
localization can be found in [21], [22]. We truncate ei 2 as —e\ = e 2 = ih, where K is a positive real 
parameter. Consequently, the formula becomes a g-series by the truncation, where q = e~ Rh . 
The fixed points in A4 n are labelled by partitions of n. The equivariant A-genus takes the 



following form at the partition A of n: 

(2mR)- 2n A T2 (Rt_ im , M n )\ x = (-) n (^y n U\ _M^ m \ > ( 2 - 26 ) 



(2mR)- 2n A T2 (Rt_ lh , lh , M n )\ x = (-) n ^y n (J[h(s)y q ^s x (qn\ (2.27) 



s x ( q n = q nW+ ^u( 1 - ( i his] T > ( 2 - 28 ) 



h \ 2n TT Mf) 

11 1 _Mii 

sex U 2 -?" 

where /i(s) denotes the hook length of the box s of the Young diagram A. The right-hand side 
of the above formula can be expressed by using a Schur function as 

ft \ 2r V-rr, , A 2 ^A_ 
sex 

where s\(q p ) is the Schur function s\(xi, Xi-, • • ■ ) specialized to x, = q l ~^ , and «(A) = 2 ^V j\ eX (j- 
i). To obtain ( 12.27ft . we have made use of the g-hook formula [23] 

sex 
where n(X) = J2i>i(i ~~ l)^i- Actually, taking account of the formula 2n(A) + |A| = — k(A)/2 + 
^ seA /i(s), the g-hook formula gives rLeA! 1 ? - ^ - — 9^~ ) = Q^ r s x(q p )- By plugging this 
into the right-hand side of (12.261) . we obtain (12.271) . 

Similarly, the fixed points in Ai n x M 4 are given as (A, 0), where A is a partition of n and 
is the origin of 1R 4 . The equivariant Chern character takes the form Tr e~ kR:F - ih ' ih I = Ok(X), 
where Cfc(A) is given by 

oo 

O k {\) = (1 - q~ k ) J2{q KXi ~ i+1) ~ q ki ~ i+1 A + 1 . (2.29) 

i=X 

The above functions have been exploited in [5], [23] from the Ad gauge theory viewpoint, with q 
or q k being replaced by a generating spectral parameter. 

By the localization formula [25], integration in the right-hand side of ( 12.241) reduces to a 
summation over contributions from the isolated fixed points of the T 2 -action onM n x (M 4 x 
• • • x M 4 ). The contribution from the fixed point consists of two factors. One is the value of the 
integrand at the fixed point, and the other is the Jacobian factor in a change of coordinates 
in a neighborhood of the fixed point. By multiplying these two factors and summing up them 
over the fixed points, we obtain the right-hand side of (12.241) . The fixed points are of the form 
(A, (0, • ■ ■ ,0)) G M n x (IR 4 x ■ • • x IR 4 ), where A is a partition of n and (0, • ■ ■ ,0) is the origin 
of 1R 4 x • • • x IR 4 . The values of the integrand at the fixed points can be expressed by using 



the formulas (12.271) and (12.291) . The Jacobian factor can be further factorized by the product 
structure of Ai n and K 4, s. The Jacobian factor originating in the change of coordinates of Ai n 
in a neighborhood of A takes the form (2ii/fij (Ylsex M s ))~ 2 ' while the equivariant volume of 
M 4 x • • • x M 4 is given by the products of H~ 2 for each Oel 4 . Therefore, the formula (12.241) 
becomes eventually the statistical sum over partitions given by 

(wi o k y h ' h =H» yi /^A(g') 2 n^ 2( ^( A )> (2-30) 

\ Jw& I n—instanton 

a A : partition of n a 

where the summation is taken over partitions of n. 

Although we have not taken into account, the Chern-Simon term can be added to a 5d 
gauge theory, with the coupling constant being quantized, in particular, for the U(l) theory, 
m — 0, ±1. It modifies the right-hand side of ( 12.301) by giving a contribution of the form 
(— ) m l A lg 2 ( for each A [TJ. Hereafter, we consider the case of the U(l) theory having the 
Chern-Simon coupling, m = 1. The corresponding generating function (12.251) becomes 



Z% h {t) = J2Q lXl s.( q n 2 e h - 2 ^^ W , (2.31) 

A 

where the summation is taken over all the partitions. 



3 Integrability of 5d J\f = 1 SYM in Q background 

We can regard the generating function (I2.3ip as a g-deformed random partition. To see this, 
note that the 4d limit R — > makes q = e~ m -^ 1, therefore, by using ( 12.281) . the Boltzmann 
weight Q w sx(q p ) 2 in fl2~3TD takes at this limit, the form (A/fr) 2|A| (]l seA Hs))' 2 , which is the 
standard weight of a random partition, that is, a Poissonized Plancherel measure of symmetric 
group. It can be also viewed as a melting crystal model, known as random plane partition. 
The corresponding model is studied in [13] as a melting crystal model with external potential, 
where the Chern characters Ok in ( 12.311) correspond precisely to the external potentials. 
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3.1 Melting crystal model with external potentials 

A plane partition tt is an array of non-negative integers 



fill VTi2 7Ti 3 
7T 2 l 7T 2 2 VT 2 3 
7T31 7T 3 2 7r 33 



(3.1) 



satisfying 7Ty > 7r^ +1 j and 7Tjj > 7ry+i for all i, j > 1. Plane partitions are identified with the 
3d Young diagrams. The 3d diagram 7r is a set of unit cubes such that TCij cubes are stacked 
vertically on each (i, j) -element of n. Diagonal slices of tc become partitions, as depicted in 




(a) 




(b) 



Figure 1: The 3d Young diagram (a) and the corresponding sequence of partitions (b). 



Fig.l. Denote 7r(m) the partition along the m-th diagonal slice, where m G Z. In particular, 
7r(0) = (vtii,7T22, • • •) is the main diagonal one. This series of partitions satisfies the condition 



-< tt(— 2) -< vr(-l) -< tt(0) >- ?r(l) >- tt(2) >- 



(3.2) 



where /i >- v means the interlace relation; ft >~ u <^=^> [i\ > V\ > \ii > v<i > jjl^ > ■ ■ ■ . 

The Hamiltonian picture emerges from the above interlace relations, by viewing a plane 
partition as evolutions of partitions by the discrete time m. In particular, transfer matrix 
formulation using 2d complex free fermions is presented in [26] , by taking advantage of the 
well-known realization of partitions in the fermion Fock space. Let ip{z) = 'Yl im£ i'4 , m z ~ rn ~ 1 and 
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^i 2 ) = Emez C z_m tie complex fermions with the anti-commutation relations, {/0 m , ^*} = 
5m+nfl and {ip m ,ipn} = {'Pm^n} = 0- Partitions can be realized as states of the fermion Fock 
space. For a partition A, the corresponding state reads 

oo 

|A)= ±H^ Xi ^% +1 \0) , (3.3) 

8=1 

where |0) denotes the Dirac sea with the vanishing U(l) charge, and is defined by the conditions 

ip m \0) = for Vm > , ^41°) = ° for Vm > 1 . (3.4) 

We may separate the interlace relations (13. 2p into two parts, each corresponding to (forward 
or backward) evolutions of partitions for m < and m > towards A = 7r(0). These two types 
of the evolutions are described in the transfer matrix formulation by using operators G± of the 
forms 

OO k 

G ± = exp{^^ I ^-^J ±fe }, (3.5) 

where J± k = J2nez '■ V'ifc-nV'n : are tlie modes of the U(l) current : ip(z)ip*(z) := J2 m& z JmZ~ m ~ l ■ 
In accord with the types of the evolutions, G± generate partitions from the Dirac sea as 

(0|G + = J>(g")(A|, G_|0) = ^ SA (/)|A). (3.6) 

A A 

In the free fermion description, we can convert the loop operators O k to fermion bilinear 
operators O k of the form 

oo 

6 k = (l-q- k ) J2 Q kn ■■ ^-n^n : +1 , (3.7) 

n=— oo 

Actually, the state (13.31) becomes the simultaneous eigenstate and reproduces the Chern char- 
acters as the eigenvalues: 

6 k \\) = O k {\)\\) . (3.8) 

Therefore, taking account of (13.61) and (13.81) . we can express the generating function (12.311) in 
the fermionic representation as 

-. oo 

Z U JS,in(t) = <0|G + Q L °exp{-5>a fc }Gqo>, (3.9) 

fe=i 

where L = J2 n ez n '• 'i'-n^n '■ ls a special element of the Virasoro algebra. 
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3.2 The integrable structure 

The fermion bilinears Ok can be regarded as a commutative sub-algebra of the quantum torus 
Lie algebra realized by the free fermions [33]. The adjoint actions of G± on the Lie algebra 
generate automorphisms of the algebra. Among them, taking advantage of the shift symmetry, 
we can eventually convert the representation (J3.9P into the expression [T3] 



1 oo 

Z-SdV = (0|exp{^^(-) fc (l-g- fe )t fe J fc } 

fe=i 

-. CO 

x g^« e xp{— £(-) fe (l-g-%J_ fe }|0>. (3.10) 

fe=i 

In the above formula, g* is the element of GL(oo) given by 

g 5dl/(l) = q f GG+Q L 0GG+q f ^ (3>11) 

where W 7 = PVg = X^nez n2 '• i'-ni'n '■ is a special element of W^ algebra. The loop operators 
Ok are converted to J/, or J_k in the formula (13.101) . Actually, J±k eventually become equivalent 
in the formula, since g* satisfies the property [13] 

J k gl dm = gl dU{l) J-k , for k > 0. (3.12) 

Viewing the coupling constants t as a series of time variables, the right-hand side of (I3.10p 
is the standard form of a tau function of 2-Toda hierarchy [27]. However, by virtue of (I3.12p . 
the two-sided time evolutions of 2-Toda hierarchy degenerate to one-sided time evolutions. 
This is precisely the reduction to 1-Toda hierarchy. Thus the generating function becomes a 
tau function of 1-Toda hierarchy. The Toda hierarchy also has a discrete variable s, which 
corresponds to the U(l) charge of the Dirac sea. The formula (I3.10P is easily generalized to the 
cases, where the charge s is interpreted as the vev of the Higgs field. 

4 Integral representation of energy functional 

A partition A = (Ai, A 2 , ■ • • ) gives rise to a series of decreasing integers \ — i {i = 1, 2, • ■ • ). 
Vice versa, a series of decreasing integers xi > x% > • ■ ■ , in which the j-th integer becomes 
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Xj = —j for j > 3i Q , gives a partition. Such a series of decreasing integers is described by its 
density function as 

oo 

Px(x) =J2S{x-{Xi-i)), xeR. (4.1) 

t=i 

The statistical sum in the right-hand of ( 12.3ip can be converted to a statistical sum over the 
density functions. 

2-T,»M = £^<>], (4.2) 

Pi-) 

where £[/?(•)] i s the energy functional obtained by taking the logarithm of the statistical weight. 

E[px{-)\ = -\og{Q^s x ( q n 2 e h - 2 ^^°^} . (4.3) 

The density function is rather implicit in the above expression of the energy functional. 
However, taking the consideration in pQ, we can convert the expression into an integration of 
a quadratic form of the density function. Prior to giving such an integral representation of the 
energy functional, let us describe the basic formula as follows: Let f(x) be a function. For a 
partition A, we put 

sex 
where h(s) denotes the hook length of the box s G A. The sum over the boxes in the right-hand 
side of (14.41) can be converted to an integration of a quadratic form of the density function. To 
see this, we take a function g(x) which satisfies 

g(x + l)-2g(x)+g(x-l) = f(x), (4.5) 

9(0) = 0. (4.6) 

By using such a function g(x), we can eventually express F\ as 

F x =\( dxdy g(\x - y\) A Px (x)A Px {y) , (4.7) 

where A denotes a difference operator of the form 

Ap(x) = p(x) - p(x - 1) . (4.8) 
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We provide a proof of the formula (14. 7p in Appendix O 

As an application of the above formula, we describe an integral representation of the loga- 
rithm of the hook polynomial 

Hx(q) = I^ 1 " Q HS) ) ■ ( 4 -9) 

Note that, by taking f(x) = log(l — q x ) in (14.41) . we find F\ = logH\(q). Therefore, by choosing 
a function g(x ; q) to satisfy the conditions ( 14.51) and (I4.6p . where f(x) = log(l— q x ), the formula 
(14. 7p gives the integral representation 

\ogH x (q) = -/ dxdy g(\x-y\; q)Ap x (x)Ap x (y). (4.10) 

In the above formula, we can take g(x ; q) to be the logarithm of the g-analogue of the Barnes 
G-function as 

e 9(x ' q) = (l-g)^G 2 (x + l;g). (4.11) 

To see that the above g(x; q) actually satisfies (14.51) and (14. 6p . we sufficiently note that G 2 (x; q) 
is the second cousin in the hierarchy of the multiple g-gamma functions G n (x; q) (n = 0, 1, • • • ), 
which are defined by the following conditions [28J: 

00 G n (x + 1; q) = G n -i(x; q)G n (x; q) , 
(ii) G n (l;q) = 1, 

rf n + l (4-12) 

(»«) j-^riogG n +i(^ + i;?) >o, x>o, 



(w) G (x;g) = [x] 



</ • 



where [x\ q = \z~- The infinite product representation of G2(x;q) follows from the above 



-a 

conditions as 



G 2 (x + l;q) = (l-?r^n{(^^)V-<? fc r}- (4-13) 

fc=i ^ ^ ' 

The above infinite product representation leads to the expansion of g(x; q) in positive powers 

of q x as 

°° 1 ( x 1 1 °° q nx 

fa; q) = 2^ - [ 1 _ q - n ~ (i_ g n)(i- g -«)| + A. n (l-q»)(l-q- n ) ' (414) 
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Let us describe an integral representation of the energy functional (j4.3p . For this end, we 
first divide £ into two parts 

£[p x {-)\ = S 1 [p x (-)]+£ 2 [p x (-)}, (4.15) 

where 

£ 1 [ Px (-)] = -logJQ'^Of) 2 }, (4.16) 

oo 

E 2 [px(-)] = -h- 2 J2t k O k (X). (4.17) 

fe=i 

To obtain an integral representation of £i, note that the g-hook formula (12.281) enables us to 

factorize Q^s x (q p ) 2 into the product q~ K ^/ 2 x ( Ylsex Ql h ^ ]H x (q)~ 2 , where the logarithm of 

the first factor reads as 

,,s ]ng- n /" + °° 

-logg-^ = — dxx 3 Ap x (x), (4.18) 

6 J-oo 

while an integral representation of the logarithm of the second factor is easily obtained from 
the previous representation of logi^^)- By summing up these two, we eventually obtain the 
integral representation of £\ as 



£M-)] = jp I dxx'Ap(x) 



-oo 



lo g<? / ,3, 



+ / dxdyg(\x-y\;q,Q)Ap(x)Ap(y), (4.19) 



where g(x; q, Q) is the logarithm of the following combination: 

! r^\ x(x— 1) x(x — 1) , v x(x — 1) , , 

e g(x;i,Q) = q r J - q -- L vr J -(l-q)-hr i G 2 ( x +l;q). (4.20) 

By using (14.131) . the expansion of g(x; q, Q) in positive powers of q x takes the form 

x(x 2 — 1) 



g(x;q,Q) = y^-i-, —~T\ n^f\ ^a f a log( 5 



12 



logg 



— nx 

+ Y" - (4:21) 

^ n(l - q n )(l - q- n ) ' K ' J 



As regards the second part, we note 



/+oo p+oo . . 

dxq kx Ap x (x) = - dxq kx (p x {x)-p x {x-l)) 

-oo J —oo 
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oo 



The right-hand side of this formula becomes a finite sum by cancellation of terms between the 
two sums and gives rise to Ok{\). Thus we find 

/+oo 
dxq kx Ap x (x). (4.23) 

-oo 

Therefore, the integral representation of £2 becomes 

/+00 °° 
dxJ2tkq kx Ap(x). (4.24) 

-°° fe=i 

5 Thermodynamic limit and Riemann-Hilbert problem 

We consider the field theory limit of the U(l) theory, which is achieved by letting h — » and 
amounts to the thermodynamic limit of the melting crystal model or the g-deformed random 
partition. Actually, the statistical average of the number of boxes of a partition , where the 
partition is identified with the Young diagram in the standard manner, becomes of order h~ 2 . 
This implies that, as h goes to zero, partitions that dominate are those of order fr 2 . To realize 
the thermodynamic limit, we have to rescale the Young diagrams by changing the size of each 
side of a box H times the original one. Correspondingly, the rescaling of partitions can be 
organized in terms of the density function, by rescaling x to u = hx as 

p(x=j:) =p®(u) + 0(h), (5.1) 

where p(°>(u) denotes the scaled density function. 

The statistical sum over the density functions in ( 14.21) can be replaced, as H goes to zero, 
with the sum over the scaled density functions 

*SU«)= E'-*"*" 11 , (5.2) 

p(0)(.) 

where £^ denotes the classical energy functional obtained from £ by 

£[p(x=^)} = ^{£ {0 V 0) (u)}+O(h)}. (5.3) 
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Taking account of the expression (j5.2p . the field theory limit is the semi-classical limit realized 
by minimizing the classical energy functional. 

The classical energy functional can be obtained from (14.191) and (I4.24p by scaling the density 
function p and the parameters as well. Noting that the difference operator A becomes the 
differential operator hd/du as fi goes to zero, the classical energy functional takes the form 

£W[ p i°)(.)] = / dudvg {0 >(\u-v\ ;R,A) , v ' H - v ; +— / duvT , 

Ju+v du dv 6 J^ du 

J-co k=1 au 

where g^\u\ R, A) denotes the classical limit of g(x; q, Q) as 

g(x=^, q = e- m ,Q = (RA) 2 ) = ±{g(°\u;R,A) + 0(h)}. (5.5) 

Explicitly, using the formula (14.211) . it is given by 

(0) R 3 logflA 2 C(2) ,C(3) 1 ^e- nRu 
g^(u;R,A) = —u —u - —u + — - — ^— —. (5.6) 



12 2 R R 2 R 2 ^ n" 



The above function is subject to the following conditions: 



d 2 9 i0) (u;R,A) . /sinhfx 

^ = log l^|H ' (5 - 7) 



2 



d g <°)(0;fi,A) = _ 

era 



(5.* 



Obtaining the thermodynamic limit reduces to solving the saddle point equation 5S^ ' /8p^> (u) 
= 0, imposing the constraints [1] originally followed by Ap\. The variational problem can be 
eventually summarized as the following integral equations: 

PP £™dv2d u g(°\\u-v\-,R,A)^^ = -!*£.- V'(u) on supp(^) , (5.9) 

du F / ; = -1, 5.10 

j du 

+ °L„« = 0, (5.11) 

du 



where the integration symbol in the first equation means principle part, and we write the 
potential as 

oo 

V(u) = Y,^ e ~ kRU - ( 5 - 12 ) 

fe=i 

Solution of the variational problem gives the classical energy E* = £^[pi (■)], where 
p; denotes the solution. The vev of the loop operators 0& can be realized by the partial 
differentiations of Si (t) with respect to £&. Actually, it is the vev at non- vanishing coupling 
constants t. We can express this quantity by using pi in the form 

afi ° ,(() - hm( 0i ) = fie— d 4<^. (5,3) 



dtk fr-*o ' J ^^ du 

To see the above formula, note that eq. (15.91) can be integrated to 

PP r°°dvg^(\u-v\;R,A)^i^- 

-oo dv 



(0) 



= ~\ (^- + V(u)) + \c(t) on supp(^) , (5.14) 

where C(t) is the integration constant. By using the above equation and eq.( l5.10l) . we can 
express bi as 



(0)a.,W„(°)/ 



Si°\t) = -[ dudvg^\\u -v\;R, A) df> * {u) df> * (v) - C(t) . (5.15) 



u+v du dv 

The partial differentiation of the right-hand side of the above equation with respect to £& leads 
to the formula (15. 13ft . Actually, making use of eqs. (15. 14ft and ( 15.101) . it can be computed as 

d tk £i°\t) = -2 jT^^Mft/pp jy vg ^\\u-v\;R,A)^^\-d tk C(t) 

°dud tk V(u)^j^-. (5.16) 

du 

Thus, taking account of (I5.12p . we obtain the formula (15.131) . 
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The variational problem can be reformulated as the Riemann-Hilbert problem to find out a 
certain analytic function. To see this, consider the integral transform 

$(z) = r^ ducot\i-{z-u) dp \ (M) zeC. (5.17) 

J-oo 2 du 

The above function is actually an analytic function on the cylinder C* = C/^pZ, since it is 
periodic with period 2m /R. The inverse-transform is realized by the imaginary part as 

«^| 9<t( „ ±i 0). ( 5 ,S) 

Therefore, taking account of the integral transformation (15.171) . the Riemann-Hilbert prob- 
lem is stated to obtain an analytic function on C* — /, where I is an interval in the real axis, 
and satisfies the following conditions: 

K$(u±zO) = -l-—V"\u) when we/, (5.19) 

R 

9f$(u±i0) = whennGl-7, (5.20) 

$(z) -► Tl as $lz -f ±oo , (5.21) 

(I) dz z&(z) =0 (C : a contour encircling / anticlockwise) . (5.22) 

Among the above conditions, the first two correspond to the saddle point equation (15.91) . More 
precisely, they are a twice differentiated form of the saddle point equation. The last two 
equations respectively amount to (j5.10p and (15.111) . 

6 Extended Seiberg-Witten geometry of 5d theory 

The solution of the forgoing Riemann-Hilbert problem was obtained in |29j,[8] in the case 
where all the coupling constants vanish. By generalizing the argument there, we can solve the 
Riemann-Hilbert problem (I5.19p - (I5.22I) . even in the case of non- vanishing coupling constants. 
Let us employ the following curve which describes a fourth punctured CPi: 

Cfs : y + y- 1 = ^(e- fe -/3), ^C, (6.1) 

where (3 is a real parameter satisfying the condition (3 > 2RA. The above curve is an analogue 
of the so-called Seiberg-Witten hyperelliptic curve. Actually, it is a double covering of the 
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Figure 2: Ca is a double cover of cylinder. 



cylinder C*, and y is uniformized on the curve. See Fig. 2. In particular, the branch points of 
y are located at z = — -^ log(/? ± 2RA). Thereby, y has a single cut along the interval / in the 
real axis, where I = -± log(J3 + 2RA) , -± log(/3 - 2RA) . 

We make an ansatz on the solution of the Riemann-Hilbert problem (I5.19p - fl5.22p as 

1 dlogy 



§(z) = -l + (2 + tp(z)) 



R dz 



(6.2) 



where tp{z) is a certain analytic function on Cp. Through the above ansatz, conditions (15.191) - 
(15.221) impose constraints on (p(z). Actually, the first three conditions are translated as 



± — V'"{u)JP(u-iO) uel, 
S(p(u±iO) = <j lR 

o ueR-i 

<f(z) = 0(1) as 3ftz — > +oo , 
<p(z) = 0(e Rz ) as^^-oo. 

In the above, P(z) is a quadratic polynomial of e Rz given by 

1 dz 

—diogy = — ; 
R ^Pjzj 

In other word, 



(6.3) 

(6.4) 
(6.5) 



(6.6) 



P{z) = (1 - (3e Rz ) 2 - {2RAe Rz ) 2 
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(6.7) 
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Figure 3: The cut I and the contour C on the Riemann sheet of C 



0- 



An analytic function satisfying the above three conditions can be realized by a contour 
integral of the form 



<p(z) = — (f dw( I + coth -(z-w))v'"( w)JW\ 
Aiii Jq \ 2 / 



w) 



(6.* 



where C is a closed curve that surrounds the interval I but leaves z outside as depicted in Fig. 
3. To see that the above satisfies the conditions (16.3p - (16.5l) . note that the contour integration 
in the right-hand side of (16.81) can be written as an integration along / of the form 



(p(z) = -4 fdv(l + coth ^(z - v))v'"(v)y/P(v - iO) . (6.9) 

2m J i \ 2 J 

When z = u ± iO, taking account of (16.71) . the imaginary part of the right-hand side of this 



equation becomes ±V'"(u)^P(u — iO)/iR if u E I, otherwise vanishes. Therefore, tp(z) given 
in ( 16.81) obeys the condition ( 16.31) . The constant term of the kernel function in the contour 
integral (I6.8P is required in order to satisfy the boundary conditions (16 .4p and (16. 5p . To see 



this, note that the kernel function has the asymptotics 



R 



1 + coth -z = 2j2 e ~ nRz as3fte-+oo 



n=0 

oo 



-2^e nRz as3?z^-oo 



(6.10) 



n=l 



By plugging the above into the right-hand side of (16.91) . we find that (p(z) satisfies the conditions 
(E3D and (1631) . 
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6.1 Determination of j3 

In order to obtain the solution, taking (p(z) as in (16. 8p . we further need to impose the condition 
( 15.221) . This turns out to be an equation which determines the parameter /3 of the curve (16. ip 
in terms of t. 

We can calculate the right-hand side of (16. 8p by residue calculus. To see this, we change 
the integration variable from w to W = e~ Rw as 

**> = ^ £ W(W-e^) &"Wy/W - fl' - i^f . ("I) 

where the contour C, which is the image of C by the mapping w i— > W, encircles the interval \j3— 
2RA, /3 + 2RA] anticlockwise. In the outside of C", the integrand has poles at W = e~ Rz , oo and 
is holomorphic elsewhere; the origin W = is not a pole because V'"(W) = YlT=i^k(kR) 3 W k 
has a zero at W = 0. We can deform C outward to small circles that surround the two poles 
clockwise. Thus the contour integral reduces to a sum of the residues (times —1) at those poles, 
so that we have 

<p{z) = ~V"'(z)VpJz)+M(z). (6.12) 

The first and second terms in the right-hand side are the contribution of the residues at W = 
e~ Rz and ^ = 00, respectively. More explicitly, 



M(z) = J2t k M k {z), (6.13) 

fc=i 

k 

M k (z) = -R 2 k 3 J2 d *-ne~ nRz , fc = l,2,-.., (6.14) 

n=0 

where d n = d n {(3) denote the coefficients of expansion of y/P(z) in positive powers of e Rz : 

00 
^/pjz) = ^d n {(3)e nRz as3?^^-oo. (6.15) 

n=0 

The first few terms read do(f3) = 1 , d\{f3) = —(3 , • • • . Let us note that this expression of M k (z) 
and M(z) readily implies that 



R 2 k*e- kRz ^f¥{z) + M k {z) = 0(e Rz ), (6.16) 

-vV'"(z)VP(z) + M(z) = 0{e Rz ) as3?z^-oo. (6.17) 

R 
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We can thus directly confirm that the right-hand side of (16.12 h satisfies the boundary condition 
(16. 5p . Actually, to solve the foregoing Riemann-Hilbert problem, we could have started from 
(16.121) and proceeded to verifying the other conditions. 

We now write down $(z), plugging the expression ( 16.121) into ( 16.21) . as follows: 

By using this formula, the condition (15.221) becomes 

<bz(2 + M(z))dlogy = 0. (6.19) 

This equation can be understood as an equation that determines the parameter (3 of the curve 
(16. ip . To write down such an equation explicitly, let us introduce the quantities 

J n {0) = ize~ mRz d\ogy, m = 0,l,2,---. (6.20) 

All the contour integrals that arise from components of M(z) in § G z(2 + M(z))dlogy, taking 
account of (I6.13P and (16.141) . have the form (16.201) and are consequently expressed in terms of 
J m {P) and d n (P). Eventually, they are arranged to provide the following equation for p: 



oo / oo 



2-R 2 J2 k\d k ((3) J J (J3) = R 2 ^ ^ kH k d k . m {(3) J J m (J3) . (6.21) 

v fc=l / m=l \k=m / 

The solution of the Riemann-Hilbert problem ( 15. 19l) - (15.22f) can be summarized as follows: 
Take the curve C / g =/ g( i ), where (3{t) is the solution of equation ( I6.2ip . Then, the analytic function 
in (16.181) gives the solution. 

6.2 The case of single coupling constant 

Let us examine the case of t — (t 1 , 0, 0, ■ • • ). In this case, eg. (16.211) is simplified to be 

(2 + R\P)J (J3) = R%JM , (6.22) 

where Jo,i(/3) are the integrals ( I6.20p . Contour integrations that appear in ( I6.20p can be eval- 
uated by computations utilizing the classical Jensen formula, and one eventually obtains 

j, os 27U , - y//3 2 - (2jAg x 

MP) = x log( ^^ ) , (6.23) 
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Figure 4: The graph of j3 = (3(tx, 0, 0, • • •). Each corresponds respectively to the cases of 
R = 0.1, R = 0.2 and R = 0.3, where A is fixed to one. 



JM 



2m {P-y/(P-(2RA)< 

-rVH 2(Mp 



13 - ^(3 2 - (2.RA) 2 



(6.24) 



We provide a derivation of the above formulas in Appendix [Dl By plugging these values of the 
integrals into (J6.22]) . eq. (16.211) in the case of t = (ti, 0, 0, • • • ) becomes 



h 



R2 P- VP 2 - ( 2RA Y 



log 



'/3 -^/p*-(2RAy 
2(i?A) 2 



(6.25) 



When the coupling constant t\ vanishes, the above equation is further simplified to be 
(3 — \/P 2 — (2RA) 2 = 2(RA) 2 . Taking account of the constraint (3 > 2RA, which is required 
to have a single cut along the real axis, one has a unique solution (3 = 1 + (RA) 2 only when 
< RA < 1. This is precisely the solution obtained in [29], [Sj. Assuming that < RA < 1 as 
well, one also has a unique solution (3(t\) of eq. (16.251) such that /3(0) = 1 + (RA) 2 . See Fig. 4. 
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6.3 Seiberg-Witten differential 

The vev of the loop operator Ok can be represented by using an analogue of the so-called Seiberg- 
Witten differential. Let $(2) be the foregoing solution described in the previous subsection. 
The inverse transform (J5.18P gives rise to the minimizer pi as 



dpf\u) 
du 



1 /„ „,/ s^dlogy(u±iO) 

=F— 2 + Mw ^ '- uel, 

2m y du (6.26) 

ueR-I. 



The integral expression (15.131) of the vev of Ok can be organized, by plugging the above formula 
into the expression and integrating by parts, to the following contour integral: 

' - lim(C fe ) = (bkRe- kRz dS, (6.27) 



dtk h-^o 27ri j q 

where dS = S'(z)dz is an analogue of the Seiberg-Witten differential. S'(z) is given by the 
indefinite integral 

S'{z) = - f (1 + ^M(z))d log y = - log y + N{z) ^/P{z) + const. , (6.28) 

where N(z) is a holomorphic function with an expansion similar to (16. 13D . Just like $(2), one 
can characterize S'(z) by a Riemann-Hilbert problem, though we omit details. 

The contour integral in the right-hand side of (I6.27P can be converted to a residue integral. 
Actually, using coordinate Z = e~ Rz , we obtain 

—y^- = \im(O k ) = Res z= JkRZ k dS) . (6.29) 

dtk n ~^° V / 

The above expression shows that Si (t) may be interpreted as a dispersionless tau function 
[30J. This is natural because the partition function is substantially a tau function of 1-Toda 
hierarchy and £; (t) gives the leading order part of the ^.-expansion of \og2_U ih (t). 

7 Conclusion and discussion 

We studied an extension of the Seiberg-Witten theory of 5d M = 1 SYM on M 4 x S 1 . We 
investigated correlation functions among the loop operators. These operators are analogues of 
the Wilson loop operators encircling the fifth- dimensional circle and give rise to the physical 
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observables of topological-twisted 5d M — 1 SYM in the Q background through the equivariant 
descent equation. The correlation functions were computed by using the localization technique. 
Generating function of the correlation functions of U(l) theory equals a statistical sum over 
partitions and reproduces the partition function of the melting crystal model with external 
potentials, where the loop operators are converted to the external potentials of the melting 
crystal model. This eventually shows that, by regarding the coupling constants of the loop 
operators as a series of time variables, the generating function is a r function of 1-Toda hierarchy. 
The 1-Toda hierarchy therefore describes a common integrable structure of 5d M = 1 SYM in 
the Q background and melting crystal model. 

The thermodynamic limit of the partition function of this model was studied by applying an 
integral formula of the energy functional, and was reformulated as a Riemann-Hilbert problem 
to obtain an analytic function on C* — /, where C* is a cylinder and I is an interval in the 
real axis, and satisfies suitable conditions. We solved the Riemann-Hilbert problem that deter- 
mines the limit shape of the main diagonal slice of random plane partitions in the presence of 
external potentials, and identified a relevant complex curve and the associated Seiberg-Witten 
differential, where the vev's of the loop operators are expressed as residue integrals of the 
differential. 

Lastly, let us briefly discuss the issue of Ad limit (see also the discussion on this issue in 
our previous work [13J). It seems quite difficult to achieve a reasonable Ad (R — > 0) limit in 
the present framework. The difficulty is rather obvious in the setup of the Riemann-Hilbert 
problem. As the definition (15.171) shows, it will not be <&(z) but rather R<&(z) that turns into 
the Ad counterpart 

, MW= rj*L*£M (r.i) 

J_ 0O z-u du 

of the integral transform of p dM in the limit as R — * 0. If we multiply both hand sides of the 
first equation (15.191) of the Riemann-Hilbert problem by R and let R — ► 0, the outcome is the 
equation 

9fM> 4d (u ± zO) = Km | (-1 - ^») =0 (u E I) (7.2) 

that has no potential term. Thus, at least in this naive prescription, the external potentials 
simply decouple in the Ad limit. (As regards the constant term —1, it is natural that it disap- 
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pears, because this term originates in a 5d Chern-Simons term [7].) One might argue that this 
difficulty can be simply avoided by rescaling tk's as tk — * R~ 2 tk before letting R — > 0. Actu- 
ally, this does not resolve the problem, because the exponential functions e~ kRu in the external 
potentials themselves become constant functions in this limit, so that the Riemann-Hilbert 
problem still takes a degenerate form: 

m 4d {u±t0) = -J^tkk 3 (uel). (7.3) 

fe=i 

Thus these naive procedures fail to generate the polynomial potentials u k , k = 1, 2, 3, ... [5] 
on the right-hand side of the Riemann-Hilbert problem. A possible remedy for this difficulty 
will be to extend the set of loop operators Ok to k < 0. One can thereby derive polynomial 
potentials from a linear combination of exponential potentials as, say, 



.,— Ru „Ru \ " 

e a \ k 



-2R 



u k (R->0). (7.4) 



This, however, causes some other problems. First of all, we have to point out that almost 
all part of our foregoing consideration assumes, implicitly or explicitly, that q is in the range 
< q < 1. Flipping the signature of k amounts to changing q — > q~ l . To introduce the 
loop operators Ok for both k > and k < simultaneously, one will be forced to consider 
the case where q = 1 (and \q\ = 1 more generally), which can only be reached as a kind of 
singular limit letting q — > 1. This is also the case for identification of the integrable structure; 
q — > 1 is a singular limit in which some properties of the quantum torus Lie algebra (crucial for 
identification of the integrable structure) break down [13]. Moreover, our method for solving 
the Riemann-Hilbert problem, too, has to be modified if the potential V(z) contains both e~ kRz 
and e kRz . Thus the issue of the 4d limit is not straightforward in our present framework, raising 
a number of interesting questions. 
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A The T 2 -action on A E 

We treat integral curves of V eie2 as curves on IR 4 in this appendix. Denote the integral curve 
passing through x by r y x . More precisely, it can be expressed as the solution of the differential 
equation 

dx(s) 



ds 



V eij£2 (x(s)) , x(0)=x. (A.l) 



Let A be a gauge potential on the S'f/(A r )-bundle E on M 4 . We consider the parallel 
transports along the integral curves. In particular, along the above 7 X , the parallel transport of 
the fibre E x i s \ to E x is described by the holonomy operator Pe Ja, M , where the path-ordered 
integration is achieved from x(s) to x along 7*. Such parallel transports provide endomorphisms 
r eit e 2 (s) of E by the formula 

(r £l , £2 (s)0)(x) =Pe-Zw A <f>(x{s)), (A.2) 

where is a section of E. The above endomorphisms define a T 2 -action on f2°(IR 4 ,.E). In 
particular, they satisfy the relation t £1j£2 (s) • 7" ei!e2 (s / ) = T ei , £2 (s + s'). The infinitesimal action 
is therefore given by 



d . , 



(A.3) 



s=0 

where t 6l)6a denotes the generator of T 2 that gives V tlt€2 . The right-hand side of (1A.3I) can be 
computed by using (1A.2J) . To see this, note that x(s) takes a form x(s) = x + sV ei:e2 (x) + 0(s 2 ) 
for a very small s, thereby we have Pe - ^*) = 1 + sV^ e2 A^{x) + 0(s 2 ) and (p(x(s)) = 
4>(x) + sV^ e d^ipc) + 0(s 2 ). By combining these two, we can write the formula (IA.2I) as 

(r £l , 62 (5)0)( ; r) = ^ + sV^ 2 A,{x) + 0{s 2 ))(^{x) + sVl e M{x) + 0{s 2 )) 

= cj){x) + s t Vtl , t2 dA(/>(x) + 0(s 2 ) , (A.4) 

where iy e e means an operation of contraction with the vector field V ei>e2 . Thus, by plugging 
the above into the right-hand side of (1A.3I) . we obtain 

t ei;€2 • (f) = L Veue2 d A (f) . (A.5) 
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Similarly, the endomorphisms induce a T 2 -action on Ae, which is the space of all the gauge 
potentials on E, by the formula 

d T ^ 2 (s)-AT eue2 {s) = T euea (s)d A , (A.6) 

where d A = d + A is a covariant differential on E. In the above formula, we use the same 
symbol to denote the T 2 -action. By using the holonomy operator, the formula can be written 
in an explicit form as 

d + (r ei , ea (s) ■ A) (x) = Pe- f *M A (d + A{x{s))) (Pe~ f: w A Y , (A.7) 

Let us describe the infinitesimal form of the above action. For a very small s, the right-hand 
side of (IA.7I) can be computed as 

-l 



1 + sV^A^x) + 0(s 2 )) {d + A(x) + s (y? uea d„Ap + A^V-^ {x)dx» + 0(s 2 ) } 



x ( I - sVl e2 A,(x) + 0(s 2 ] 



= d + A{x) + s(vi t2 d u A, - V^ e2 d,A u + [Vl e2 A u , A M ])(x)dx" + 0(s 2 ) 
= d + A{x) + s L Vei>e2 F A (x) + 0(s 2 ) , (A.i 

where F A = dA + A A A is the curvature two form. Thus, the formula (1A.7J) reads as 



r ei , ea (s) • A = A + si Veue2 F A + 0(s 2 ) . (A.9) 



Therefore the infinitesimal form of the T 2 -action becomes 

d 
1-61,62 ' A = — r £1 6 2 {s)A 

as 



s=0 

^v eue2 F A . (A.10) 



B Proof of the identity (EIS) 



We first rewrite the Q-transformations (12.61) and ( 12.71) in the forms 

Q eue2 A{t) = ij(t), 
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dA(t) 

Q*i,*M) = d A{t) (f)(t) - L Vn ^F A{t) — , 

Q eilfa <Kt) = -tv^m- (B- 1 ) 

By using the above expression, we can easily see that the combination i^Wj — ip(t)+<f>(t) satisfies 
the identity 

dA{t) 



(d A (t) ~ V n , 62 + Qe 1: e 2 ) (^(t) " l/>{t) + 0(t)) 



dt 



B.2 



The above identity may be interpreted as a loop space analogue of the equivariant Bianchi 
identity [25J. Actually, when A(t),i/j(t) and <j)(t) are constant loops, that is, when they do not 
depend on t, ( 1B.2[) reduces to 

(d A ~ L Vei , a + Qei,e 2 ) (Fa ~ j> + <P) = , (B.3) 

where the combination Fa(x) — ip{x) + <j>(x) is naturally identified with the T 2 -equivariant 
curvature of the universal connection [3j, and (1B.3J) substantially describes the equivariant 
Bianchi identity. 

Let us derive the formula (12.181) by using the identity (1B.2J) . To do this, note that, since 
ip(t) is a Grassmann-odd one-form on M 4 , while F^it) an d (f>(t) are Grassmann-even two- and 
zero-forms, we have 

(d eit62 + Q ei , e2 )W(x;t 1 ,t 2 ) 

ds W{x; h, s) A (d euea + Q ei>ea ) (f a{s) - V>(s) + 0(s)) (z) A W(x; s, t 2 ) , (B.4) 



't 2 
where W is another generalization of the path-ordered integral (I2.10p , and is given by 

't 2 
Since 0(x) is the trace of W(x; R, 0), using (IB.4I) . we can compute (rf 6l)62 + Q £lj(E2 )(9(a;) as 



W{x; UM) = Pexpj- / ' rft(F A(i) +V(t) + 0(t))(z)} . (B.5) 



(d euea + Q euea )0(x) 

dtTr{w(x; R, t) A (d ei , ea + Q 6lj£2 ) (> A(f) - </>(t) + 0(t)) (x) 

AW(a?;t,0)}. (B.6) 
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By using (IB. 21) . The right-hand side of (1B.6J) can be further evaluated as 



(d eue2 + Q eue2 )0(x) 

■R 



dtTr\W(x;R,t) A 



dA(x,t) 



dt 



AW(x;t,0)\ (B.7) 

+ / dtTrlw{x;R,t) A A(x,t) A (F A(t) -^(t) + <j>(t)\(x) AW{x;t,0)\ (B.8) 
W{x; R, t) A (,F A{t) + ip{t) + <f>(t)\ (x) A A(x, t) A W{x; t, 0) j . (B.9) 



R 



dt Tr 



In the above, we can arrange the first term of the right-hand side as 

dEZD = - 



dt Tr | — p^ A W{x- 1, 0) A W{x; R, t) 
dtTil dA ^ X,t) AW(x;t,t-R) 



o 



dt 



Similarly, the second term and the third term can be arranged into 

■R 



dESD = / dtTr{A(x,t)A(F A{t) -i){t) + (f>(t)yx)AW(x;t,t-R)}, 
(TB~9|) = - dtTr(A(x,t) AW(x;t,t-R) A (F A{t ^ R) - i/j(t - R) + <j)(t - R) 



'o 
These two are particularly combined to give 

f-R 



i •*- , 



(1B~8]) + ([13]) = - / dt Tr i A(x, t) A 



cft^x; t, t - R) 



dt 



Therefore, by summing up (1B.7|) . (1B.8|) and (IB. 91) . we find 



(d ei , e2 + Q ei , e2 )0( 



R dtTil dA ^ t) AW(x;t,t-R) 



("i 



o 



rft — 



tft 



Ir %i)AW(a;;M--R) 



(it Tr <^ A(x, t) A K -^j - 



(B.10) 



The last integral is integrated to be Tr A(x, 0) A W(x; 0, -R) - Tr A(x, R) A W(x; R, 0), and 
becomes, by virtue of the periodicity, zero. Thus we obtain the formula (I2.18p . 
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C Proof of the formula (14.71) 



The formula (14.71) is stated as follows. 

Formula 1 Let f(x) be a function on R. Choose a function g(x) to satisfy the conditions jjj4-5\ ) 
and lfi4-6\ )- Then, for any partition A = (Ai, A2, • • • ), the following identity holds: 

J2f(h(s))= f dxdyg(x-y)Ap x (x)Ap x (y), (C.l) 

sg\ Jx> y 

where h(s) denotes the hook length of the box s G A ; and Ap x (x) = p\(x) — p\(x — 1) is a 
difference of the density function (4A) of X. 



To prove the above formula, note that the density functions of a partition A and its conjugate 
partition A = (Ai,A 2 ,---)> where the Young diagram A is obtained by flipping the Young 
diagram A over its main diagonal, satisfy the relation Ap x (x) = Ap x (— x). We can thereby 
rewrite the right-hand side of ( 1C.1I) as 

/ dxdyg(x-y)Ap x (x)Ap x (y)= dxdy g{x - y) Ap x (x)Ap x (-y) (C.2) 

J x>y J x>y 

We compute the right-hand side of (1C.2I) . For this end, it is convenient to introduce the 
truth function \ by 

I 1 iff n > 1 , 

X(n) = { ~ (C.3) 

j otherwise . 

Owing to the ^-functions in the density functions, the right-hand side of ( 1C.2I) can be expressed 
by using the above x as 



dxdyg(x - y) Ap x (x)Ap x (y) 



x>y 
00 



J2 gi^i + Xj-i- j)x(Ai + Xj - i - j) + J2 9(^ + Xj-i-j + 2) X (Xi + Xj-i-j + 2) 

i,j=l i,j=l 

00 

-2 Y, 9(\ + A,— i - j + l)x(\ + ~X j -i-j + l). (C.4) 



We can further simplify (IC4j) . To see this, note that the combination Aj + Aj — i— j becomes 



< —2 when the pair (i,j) is not a box of the Young diagram A, while it becomes > when the 
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pair is a box of A. More precisely, when the pair (i,j) is a corner of A, which means that (i,j) 
is a box of A but neither (i + l,j) nor (i, j + 1) is, the combination Aj + Xj — i — j is equal to 
0, otherwise it takes > 2. Therefore, (IC.4I) is eventually simplified as 

dxdyg(x - y) Ap x (x)Ap x (y) 

x>y 

= J2 {^ + A, - ^ - J + 2) - 2g(Xi + Aj - i - j + 1) + g{\ + Xj-i- j)} 

corners of A 

In the above, by virtue of the condition (I4.5p . the first term of the right-hand side becomes 
J2(in -)e\ f(h(h j))i while the second term vanishes by the condition (14.61) . Thus we obtain the 
formula (RTT]» . 



D Derivation of eqs.(E231) and (E23) 



We derive eqs. f!6.23|) and (j6.24p by using the classical Jensen formula and its variant. Let a be 
a complex number such that \a\ ^ 0,1. The classical Jensen formula is an integration formula 
of the form [31] 

i r 27T 

— d6\og\e ld -a\ = logmax(l,|a|). (D.l) 

We also use the following variant of the above formula: 

_L f dd ( e ™ + e -*») l g | e * _ a | = J fc (D.2) 

2vr / n 1 

--^a k iff < lal < 1 
V. k 

where k — 1, 2, ■ • ■ . 

Let us first express J m (/3) (16.201) in terms of a contour integral on the ?/-plane. To do this, 
note that, taking account of the relation e~ Rz = RA(y + y^ 1 ) + /3, the contour C in the right- 
hand side of (I6.20p can be chosen so that it maps to the unit circle \y\ = 1 on the y-plane. We 
also note that, on the unit circle, by virtue of the condition (3 > 2RA, RA(y + y~ l ) + (3 takes 
positive real numbers. Thereby, the contour integral in the right-hand side of (16.201) can be 
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written as 



Jm((3) = -i / ^-{RA{y + y- 1 ) + p}™ log \RA(y + y~ l ) + (3\ . (D.3) 

R J\ y \=i y L J 

The above integral can be converted to a combination of the phase integrals that appear in 
the formulas ( ID. II) and ( ID. 21) . To see this, note that RA(y + y~ x ) + (3 is factorized into 



i?A(y + y' 1 ) +P = RAy-\y - a){y - a" 1 ) , (D.4) 

where a = —(3/2RA — y/([3/2RA) 2 — 1. By plugging the above into the right-hand side of 
(1D.3|) , the integral becomes a sum of phase integrals of the form 



i log RA 



R jo 

»2tt 



2- 



Jm(P) = 15— d6\RA(e M + e- w )+(3 



- dO\RA{e ie + e~ ie ) + (3\ \og\e ie - a\ 

- / d6lRA(e ie + e- w ) + (3\ log \e ie - a^\ . (D.5) 



The above phase integrals can be computed by applying the formulas ( ID. II) and (ID. 21) . For the 
cases of m — 0, 1, we particularly obtain eqs. (j6.23p and (16.24)) . 
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